Inseparability criteria based on matrices of moments 
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Inseparability criteria for continuous and discrete bipartite quantum states based on moments of 
annihilation and creation operators are studied by developing the idea of Shchukin-Vogel criterion 
[Phys. Rev. Lett. 95, 230502 (2005)]. If a state is separable, then the corresponding matrix of 
moments is separable too. Thus, we derive generalized criteria, based on the separability properties 
of the matrix of moments, are thus derived. In particular, a new criterion based on realignment 
of moments in the matrix is proposed as an analogue of the standard realignment criterion for 
density matrices. Other inseparability inequalities are obtained by applying positive maps to the 
matrix of moments. Usefulness of the Shchukin-Vogel criterion to describe bipartite-entanglement of 
more than two modes is demonstrated: We obtain some previously known three-mode inseparability 
criteria originally derived from the Cauchy-Schwarz inequality, and we introduce new ones. 



I. INTRODUCTION 

In recent years, the study of continuous-variable (CV) 
systems from the point of view of quantum information 
has attracted much interest, stimulated by experimental 
progress (see [3,0 an d references therein). In particular, 
the theory of quantum entanglement for CV systems has 
been considerably developed, including the derivation by 
Shchukin and Vogel [3( of a powerful inseparability cri- 
terion of bipartite harmonic quantum states based on 
Partial Transposition (PT) [1, H[, the so-called Positive 
Partial Transposition (PPT) criterion. The PPT crite- 
rion says that a separable state remains positive under 
partial transposition, therefore a Non-positive-Partial- 
Transposition (NPT) state must be entangled. Shchukin 
and Vogel have demonstrated that their criterion in- 
cludes, as special cases, other well-known criteria of en- 
tanglement in two-mode CV systems, including those de- 
rived by Simon [6|, Duan et al. J3, Mancini [8(, Raymer 
et al. [9|, Agarwal and Biswas [lOj, Hillery and Zubairy 
[TTI | . Thus, the Shchukin-Vogel (SV) criterion can be 
considered a breakthrough result, which shows a com- 
mon basis of many inseparability criteria for continuous 
variables (in particular, the results of Duan et al. [Tj 
seemed previously to be entirely independent of partial 
transposition). Another advantage of the SV criterion 
should be noted: it is given in terms of creation-operator 
and annihilation-operator moments, which are measur- 
able in standard homodyne correlation experiments [l2j 
(for recent reviews on entanglement detection see Refs. 

Despite the evident progresses (see also [TH, H EH, [HI, 
UtI [la ] and references therein), the theory of quantum 
entanglement for CV systems can be considered less de- 
veloped than the theory for discrete, finite-dimensional 
systems [li|. In the latter case, powerful inseparability 
criteria based on positive maps (see 1131. 19 1 and refer- 
ences therein) and linear contractions [2C 
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permutations of the indices of density matrix [24[ ) have 
been studied as generalizations of the standard PPT cri- 
terion [1, [j| . Inspired by these tools available to study 
discrete- variable entanglement, we propose a generaliza- 
tion of the Shchukin-Vogel CV approach. 

Shchukin and Vogel [|j] recognized a deep link between 
the property of positivity under the operation of PT of 
a two-mode density operator p, and the positivity un- 
der PT of the corresponding matrix of moments. In the 
present work, we obtain a more general relationship be- 
tween the separability properties of the density operator 
and of the matrix of moments. Namely, we show that if 
a state is separable, then a suitably designed matrix of 
moments is separable too. This will allow us to apply 
all known separability criteria (not only the PPT one) to 
the matrix of moments rather than directly to the density 
matrix. For the sake of clarity, we will analyze explicitly 
mainly the bipartite two-mode case; anyway, the results 
can be extended to the multimode (see Sect. VII) and 
multipartite case. 

As the objectives of the paper are of wide range, let 
us first specify the main goal and results of the paper. 
We analyze the Shchukin-Vogel inseparability criterion 
for matrices of moments from a new perspective useful 
for generalizations along the lines of the standard insep- 
arability criteria for density matrices. More specifically, 
we emphasize the fact that separability is preserved by 
the mapping from states to matrices of moments. This 
more general approach leads us to propose new entangle- 
ment criteria based on realignment and positive maps, 
which lead to new inequalities directly applicable in ex- 
perimental tests of entanglement. 

In particular, in Sect. II, we present a general idea of 
separability criteria based on matrices of moments. In 
Sect. Ill, we review the Shchukin-Vogel criterion. In 
Sects. IV and V, we present our generalizations of the 
SV criterion based on the separability properties of the 
matrix of moments of creation and annihilation opera- 
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tors by referring to contraction maps (e.g., realignment) 
and positive maps (e.g., those of Kossakowski, Choi and 
Breuer). A few examples illustrating the applicability of 
the new criteria are shown. In Sect. VI, we discuss de- 
tection of entanglement by expressing the entries of the 
density matrix in terms of the moments. In Sect. VII, we 
briefly discuss the use of the criteria to analyze bipartite- 
entanglement of more than two modes. Finally, we give 
our conclusions. 



II. SEPARABILITY OF STATES AND 
MATRICES OF MOMENTS 

Consider two modes A and B with associated annihi- 
lation and creation operators a and a' for A and b and 

for B. Shchukin and Vogel showed that a Hermitian 
operator X = X AB is nonnegative if and only if for any 
operator / = f AB whose normally-ordered form exists, 
i.e., 

+00 

/= E c klk2hl J kl a k *tf h b h , (1) 

ki,k2,h,l2— 

it holds (/t/) x =Tr{/t/Jf}>0. 
Let us consider the operators 

fi = f£f?, (2) 

with f A = a) kl a k2 and f B = tf h b h . Here i is the 
unique natural number associated with a double multi- 
index (k, 1), with k = (fci,fc 2 ), 1 = (hjfa)- Furthermore, 
the multi-indices k and 1 are associated with unique nat- 
ural numbers k <-> {k\,k2) and I <-> (/i,Za)< Any op- 
erator / whose normally form exists can thus be writ- 
ten as / = 53j c i/i' ^ we further define the matrix 
M(X) = [Mij(X)], whose elements are given by 

M ti (X) = (f}fj) x = Tr {fjfjX}, (3) 

we have 

Lemma 1 An operator X is positive semidefinite (X > 
0) if and only if M(X) is positive semidefinite Q. 



p we have: 

Mij(p A ® p B ) 
= Tt{flf jP A ®p B } 

= Tr { (a tfel a k2 )\a^ a k ' 2 ){b ]h b l2 )^ (b^ b l ' 2 )p A <g) p B } 
= Tr { (a tfel a fe2 ) f (a^ a k ' 2 )p A }Tr {(tf h b h )\b^ b l ' 2 )p B } 
= Tr{(f A yf A p A }Tr{(f B yf B p B } 
= M A ,(p A )M B (p B ), 

(4) 

where M A k , (p A ) = Tr {(/^)t f A p A } , so that M A (p A ) = 
[M A k ,(p A )] is the matrix of moments of subsystem A in 
state p A (and similarly for B). 

A matrix of moments uniquely defines a state, i.e. if 
M(p) = M[a) then p — a. This is immediately proven by 
considering that if M(p) = M{a) then Tr {{p- o-)/t/} = 
for all fs. 

We introduce explicitly formal (infinite) bases [42j 
\k) = |k) and |Z) = |1), in which we express the matrices 
of moments: 

M( P ) - E M kW i'{p)\k){k'\ ® (5) 

kk'W 

Taking into account the one-to-one correspondence be- 
tween matrices of moments and states and (j3|, we con- 
clude 

Proposition 1 A state is separable, p = J2iPiPt ® 
pf 1 Pi > 0, YliPi = 1; */ an d on M) if t ae cor- 
responding matrix of moments is also separable, i.e., 
M (P) = Y.^ MA {pt) ® M B (P?) with M A (p A ) = 
Efcfc' M A k ,{p A )\k)(k'\ and analogously for M B (p B ). 

Notice that the local matrices of moments 
M A ^ B \pf^ B ^) in the Proposition are physical, i.e., 
can consistently be interpreted as related to a local 
state. Thus, one has to take into account the subtle 
point that a matrix of moments could be separable in 
terms of generic positive matrices, but not in terms of 
physical local matrices of moments. Such a point does 
not arise when studying the entanglement of a density 
matrix: in that case, any convex decomposition in tensor 
products of positive matrices is automatically a good 
physical separable decomposition. Therefore, it might 
be that no method based on the study of separability 
properties of matrix of moments, can distinguish all 
entangled states. 



Indeed, X is positive semidefinite if and only if (f'f)x > 
for all f = J2i dfh i-e-, if and only if ^ c*c 3 M, 3 (X) > 
for all possible (ci)i = (ci, C2, . . .). In turn, this implies 
that X > if and only if M(X) = [My(X)] is a positive 
semidefinite (infinite) matrix. We will refer to correlation 
matrices as M(X) as to the matrices of moments. 

For any density operator p AB , from Lemma Q] we have 
that the corresponding matrix of moments M(p AB ) is 
positive semidefinite. For a factorized state p AB — p A <g) 



III. PARTIAL TRANSPOSITION AND 
SHCHUKIN- VOGEL CRITERION 

Let us now recall the Shchukin- Vogel reasoning 0. 
Let us first define the operation of partial transposition. 
Given a density operator 

P = E PM,k'v\kl){k'l'\ (6) 

k,l,k',l' 
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in some fixed basis (say in Fock basis), where pkik'l' = 
{kl\p\k' V) , its partial transposition (with respect to sub- 
system B) is 

P r = Pki,k'l'\kl')(k'l\. (7) 

k,l,k',V 

Partial transposition is a positive but not completely pos- 
itive [431 ] linear map which is well defined also in an 
infinite-dimensional setting. Positivity of p T is a nec- 
essary condition for separability of p [1, |j| . We rederive 
explicitly the relation between the matrix of moments of 
p and the one of the partially-transposed state p T : 

M u ,k'v{p T ) 
= Tr {(at fel a k2 )t (a^ a k * ) (6+' 1 b l *)Hb^ b 1 '* )p F } 

= Tr{(a t ' £l a' £2 ) t (a tfe ia^)((6 t/l 6 /2 ) t (6 tr i^)) T p} 
= Tr{(o t * 1 a fe, )t(o t * i o* 4 )(6 t 'i6 , »)t(6t'i6'a) / ,} 
= M k v,k'i{p), 

following from the property b T = b^ . Therefore, the ma- 
trix of moments of the partially-transposed state corre- 
sponds to the partial transpositions of the matrix of mo- 
ments of the state. Moreover, considering Lemma [TJ we 
have: 

Criterion 1 (Shchukin-Vogel [3]) A bipartite quan- 
tum state p is NPT if and only if M{p T ) = (M(p)) r 
is NPT. 

Considering the remarks following Proposition [JJ it is 
noteworthy that analyzing the partial transposition of 
the matrix of moments we are able to conclude about 
the PPT/NPT property of the states. In particular, this 
means that the only possible entangled states, for which 
the analysis of the separability properties of the corre- 
sponding matrix of moments is not enough to reveal their 
entanglement, are PPT bound entangled states [25ll26l|. 

Given Criterion [TJ there is still the problem of analyz- 
ing the positivity of (M(p)) r . Since the matrix of mo- 
ments is infinite, one necessarily focuses on submatriccs. 
Let us define Mjy(p r ) to be the submatrix corresponding 
to the first N rows and columns of M(p T ). According to 
the original work by Shchukin and Vogel 0, a bipartite 
quantum state would be NPT if and only if there exists 
an N such that det Mm{p t ) < 0. As shown in [27fl . this is 
not correct, since the sign of all leading principal minors, 
i.e., of det M/v(p r ), for all N > 1, does not characterize 
completely the (semi)positivity of matrices of moments 
which are singular. For any (possibly infinite) matrix 
M, let M r , r = (r*i, . . . ,r N ) denote the N x N princi- 
pal submatrix which is obtained by deleting all rows and 
columns except the ones labelled by n , . . . , r/v . By ap- 
plying Sylvester's criterion (see, e.g., [28|]) we find [27j : 

Criterion 2 A bipartite state p is NPT if and 
only if there exists a negative principal minor, i.e., 



det(M(p r )) r < for some r = (r\, . . . , rjy) with 1 < 
7"! < r 2 < ■ ■ ■ < r N . 

Focusing on the principal submatrix (M(p)) r , is equiv- 
alent to considering a matrix given by moments My (p) = 
Tr{///jp} only for some specific operators fi. In turn, 
this amounts to study positivity of p (or p r , when we 
consider (M(p r )) T ) only with respect to a subclass of 
operators p f (see the proof of Lemma [TJ, i.e., with 
/ = X)i=i °n /n • Hereafter, if not otherwise speci- 
fied, we slightly abuse notation and denote by / = 
(/ru /r 2 i — ! fr N ) & subclass of the class of operators ©. 
Let Mf(p r ) = (M(/)) r with / = (/ ri ,/ r2 ,...,/ rK ) 
denote the principal submatrix corresponding to r = 
(ri, r2, . . . , r/v). Criterion [5] can then equivalcntly be 
rewritten as: 

Criterion 3 A bipartite state p is NPT if and only if 
there exists f such that det Mf(p r ) is negative. 

More compactly: 

p is PPT & V/ : detM/(p r ) > 0, 
pisNPT^B/: det M f {p r ) < 0. (9) 

Notice that in general Mf{p r ) ^ (Mf(p)) r , i.e., the op- 
eration of partial transposition and the choice of a princi- 
pal submatrix do not commute. The criterion requires to 
consider submatrices of the partially-transposed matrix 
of moments, i.e., Mf(p r ), not to take submatrices of the 
matrix of moments and study their partial transposition. 
Nonetheless, also considering the partial transposition of 
a submatrix of the matrix of moments is a test for sepa- 
rability, if the submatrix is chosen in the right way (see 
Section ITVl in particular Eq. ([TT])) 

On the other hand, for any / (i.e., for any r), the mo- 
ments which constitute the entries of Mf(p r ) and Mf(p), 
when both expressed with respect to p, are simply related 
by Hermitian conjugation of the mode b. 

IV. NEW INSEPARABILITY CRITERIA VIA 
REORDERING OF MATRICES OF MOMENTS 

In this Section, we will be interested in studying the 
separability properties of the matrix of moments through 
a reordering of its elements. Indeed, apart from par- 
tial transposition, there are other entanglement criteria 
based on such reorderings. In the bipartite setting, the 
only non-trivial one which is also independent of partial 
transposition is realignment. For a state p as in ([Jj]), the 
realigned state reads 

P R = J2 Pki,k>i>\kk')(ll'\. (10) 

k,l,k',l' 

In a finite-dimensional setting, necessary conditions for 
separability can be formulated as ||p r || < 1 [U and 
\\p R \\ < 1 @, El, where \\A\\ = TrjVIU} is the 
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trace norm of A. The converse statements, ||p r || > 1 
and > 1, are therefore sufficient conditions for the 

state to be entangled. It is worth noting that ||p r || < 1, 
contrary to the realignment criterion, is also a sufficient 
condition for separability for 2x2 and 2x3 systems |5j . 

We have seen how the partial transposition of the ma- 
trix of moments corresponds to the matrix of moments 
of the partially-transposed state, leading to the SV cri- 
terion. It is immediate to define a realigned matrix of 
moments following (TTUJ). Unfortunately, there is no sim- 
ple relation between the realigned matrix of moments 
and the realigned state. More importantly, partial trans- 
position and realignment, while both corresponding to 
a reordering of the elements of a matrix, appear to be 
on a different footing as regards their applicability in an 
infinite-dimensional setting. Indeed, the partial transpo- 
sition criterion can be stated as a condition on positivity 
of the partially-transposed state/matrix of moments, be- 
sides a condition on the corresponding trace norm. On 
the other hand, the realignment condition can be ex- 
pressed only in the latter way, so that it is not suited 
to study the separability properties of a non-normalized 
(and non-normalizable) infinite matrix, e.g in the case of 
the matrix of moments. To circumvent such an issue, in 
the following we will analyze separability properties of 
properly truncated matrix of moments, opening the pos- 
sibility to deploy the power of the techniques developed 
for finite-dimensional systems. We remark that such a 
"truncation approach" could also be applied directly to 
CV density matrices, as it was done, for example, in [17| . 
but in this work we focus on the matrices of moments. 
One of the main reasons is that, as already remarked 
about SV criterion, moments are measurable in standard 
homodyne correlation experiments. 

In the SV approach, one typically refers directly to 
the total infinite matrix of moments M(p r ) (see Crite- 
rion [1} , studying positivity of its principal minors (see 
Criterion [2|) . Instead, we propose to first truncate the 
matrix of moments M(p), and then analyze with differ- 
ent criteria the separability of the truncated matrix of 
moments. Indeed, truncation is equivalent to focusing on 
(some) submatrix. The submatrix must be chosen cor- 
rectly, avoiding the introduction of artifact entanglement 
by the truncation. The truncated matrix is positive and, 
once normalized, can be considered a legitimate state of 
an effective bi- or multi-partite finite-dimensional system. 
Explicitly, consider subsets of indices 



Ia = {fc (1) 



and the corresponding projectors P A = J2kei A I^X^I 
and Pb — J2iei B 10 ('I- Then we can define a finite- 
dimensional matrix 



M 



I All 



.(p) = (P A ®P B )M{p)(P A ®P B ) (11) 



and we have that Mj A j B (p)/Ti {Mj A i B (p)} is a well- 
defined state (positive and with trace equal to one) for a 



cLa <8> ds system, which is separable if the starting state 
p is separable. Indeed, according to Proposition [TJ if p is 
separable then M(p) is separable too; moreover, a further 
local projection cannot induce the creation of entangle- 
ment. 

As we noted at the end of Section Hill any choice of 
a principal submatrix can be described as considering 
a specific class / of operators, i.e., a restricted set of 
products of annihilation and creation operators in normal 
order. Now, we are interested in the classes of operators 
corresponding to the choice of I a and 1b- This means 
we will always consider only tensor product classes of 
operators: 



/ = r®r 



7, i n K 2 



(g) 6' 2 , . . . , 6 t( i 6' 2 ) 

/ tfe (1) fc (1) 7t; (1) W (1) ^ 



(12) 



With the help of this notation, a truncated matrix of 
moments will be denoted in the following as 



M f (p) 



J2 M kltkn .(p)\kl){k'l' 



(13) 



fc.fc'e /a 
i,i'ei B 



for an operator class /, which is given by a tensor product 
of classes (as marked by tilde) . 

Elements of matrix (fT"3|) can be reordered to get entan- 
glement criteria in full analogy to those based on reorder- 
ing of the density matrix elements. Thus, we formally 
apply to Mj(p) the "partial transposition" 

[M f -(p)) r = ]T M kim ,{p)\k'l){kl'\, (14) 

k,l,k',l' 

and the "realignment" 

(M f ~(p)) R = ]T M Hmi {p)\kk'){ll% (15) 

k,l,k'M 

in complete analogy to ([7]) and (fTU)) . Let us define the 
normalized trace norms 



» vf (p) - 



\w f <p)r\\ 

V f [P) ~ Tr{M f -(p)} 



Tr{M 7 (p)} 



(16) 



It is worth noting that, because of the tensor product 
structure of /, we have 



(M 7 (p)) r = M f -(p r ) 



(17) 



for all / and all p. 

The SV criterion can now be equivalently formulated 

as 

Criterion 4 A bipartite state p is NPT if and only if 
there exists a tensor product class f , given by lilty) , such 
that Mj(p v ) is not positive or, equivalently, i^{p) > 1. 
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The Rudolph-Chen- Wu [20|, |2l| realignment criterion for 
density matrices, can be generalized straightforwardly for 
the matrices of moments as follows: 

Criterion 5 A bipartite quantum state p is inseparable 
if there exists f, such that (MAp)) R has trace norm 
||(Mjr(|C?)) R || greater than Tr {M.-(p)}. 



More compactly: 

p is separable 
p is inseparable 



yf(p) < i, 



v/ 

3/: vf{p)>l 



(18) 



In principle, the criterion (|18[) based on the realignment 
of the matrix of moments is inequivalent to the SV crite- 
rion based on PT, similarly as, for finite-dimensional den- 
sity matrices, the Peres-Horodecki criterion is not equiv- 
alent to the Rudolph- Chen- Wu criterion. 



in the standard SV approach, Mf(p r ) with / = (l,ao). 
Then one gets 



M f (p r ): 



1 



(a&t) 
(at b) (N a N b ) 



(21) 



from which the Hillery-Zubairy criterion of entanglement 
follows [111 ]: 



det M f (p r ) = (N a N b )-\(ab^)\ 2 < 0. 



(22) 



For our state, one gets M f (p r ) = [1, -1/2; -1/2, 0], 
which results in detM/(p r ) = —1/4. 

Example 2. The realignment-based and PT-based 
criteria can also detect the entanglement of partially- 
entangled states. To show this, let us analyze the state 
|^) = ^=(|00) + |01) + 1 10) ) for which negativity is equal 

to 2/3. By choosing / the same as in Example 1, one 
gets 



A. Exemplary applications of partial transposition 
and realignment 

Let us give a few examples of application of the insepa- 
rability criteria based on PT and realignment of matrices 
of moments. We recall that (Mj(p)) T — Mj(p T ) for a 

tensor-product /. 

Example 1. To detect the entanglement of the sin- 
glet state \ip) = -i=(|01) — 1 10) ) , one can choose / = 
(1, a) ® (1, b) = (1, a, b, ab) yielding the following matrix 
of moments M;(p) = [My] = [<#£>]: 



M f (p) 



1 

(at) 
<&t) 



6t) 



(a) 
(a&t) 



(&) 



(at&) 
(N b ) 
(a<N b ) 



(ab) 
(N a b) 
(aN b ) 
(N a N b ) 



(19) 



where p — \tp)(tp\, and N a — at a, N b = &t& are the 
number operators. The only nonzero terms of (|19p for 
the singlet state are: Mn — 1, M 2 2 = M 33 = — M 2 3 = 
— M32 = 1/2. Elements of [My] can be reordered, ac- 
cording to HU) and (15]), to get (M f ~(p)) r and (M f ~(p)) R 
equal to 



Mn 


M 2 i 


M13 


M 23 " 




-M n 


Mi 2 


M21 


M22 " 


Mia 


M 22 


M 14 


M 24 




M13 


M14 


M23 


M24 


M31 


M 4 i 


M 33 


M 43 




M31 


M32 


M41 


M42 


M 32 


M 42 


M 34 


M 44 . 




.M33 


M34 


M43 


M44 . 



(20) 

respectively. Thus, for the singlet state one gets the trace 
norms, defined by (116|) . greater than 1, i.e., v T ~ = v R = 

(1 + \/2)/2, as well as negative detMj(p r ) = —1/16 and 
mineigM / ~( / 9 r ) = (1 - V2)/2. It is seen that both the PT 
and realignment based criteria detect the entanglement 
of the singlet state. It is worth noting that one could 
analyze just the submatrix of the first matrix of (120|) 
corresponding to r = (1, 4). This amounts to considering, 



M f (p) 



-3 110- 
1110 
1110 




(23) 



which implies v v - — v R — 1.1891 > 1 (as well as 

det Mf(p T ) = -1/81 < 0). Thus, the entanglement of 
the state can be detected by both criteria. As in Exam- 
ple 1, we can use the submatrix of moments Mf(p T ) = 
[1, 1/3; 1/3,0], given by ([2T]) (or, which is the same, the 
submatrix (Mj(p T )) r of the partially-transposed Mj(p) 
of ([2"3"| . for r = (1,4)), which also has negative determi- 
nant (equal to —1/9) and minimum eigenvalue, given by 
(3- VT3)/6 re -0.1. 

Example 3. The realignment-based criterion is sensi- 
tive also for some infinite-dimensional entangled states, 
as can be shown on the example of superpositions of co- 
herent states, referred to as the two-mode Schrodinger 
cat states, 

W) = Af'(\ a ,-(3)-\-a,f3)), 
W) = Af"(\a,(3)- I -«,-/?)), 

which are normalized by functions A/ 7 and A/ 7 ' of the 
complex amplitudes a and j3. As actually shown in [3[, 
the entanglement of \rp") (but also of can be de- 

tected by the standard SV criterion for / = (1,6, ab), 
for which one gets a negative determinant det Mf(p T ). 
The realignment-based criterion applied to the factorized 
/ = (1, a) ® (1, &) is also sensitive enough to detect entan- 
glement of both states \tp') and \ip"). E.g., for both states 
with a = 0.3 and (3 = 0.2, one gets the trace norms for 
realignment and PT greater than one, i.e., v R — 1.1666 

and vj = 1.1783. Note again that by analyzing determi- 
nant or minimum eigenvalue of submatrix (M^(p r )) r for 
r = (1,4), given by ([2"T]) , one can detect entanglement of 
the state by handling less moments. 
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V. POSITIVE MAPS ACTING ON MATRICES 
OF MOMENTS 

In this section we generalize the SV criterion b y ap ply- 
ing the theory of positive maps (see reviews (13. Il9|). 

The standard criterion of separability for states which 
is based on positive maps says the following [U, 0] : a bi- 
partite state p is separable if and only if every positive 
linear map A acting partially (say on the second subsys- 
tem only) transforms p into a new matrix with nonnega- 
tive spectrum, i.e., 



(icU®A B )[p A "] >0 



(24) 



(For brevity, the system-identifying superscripts are usu- 
ally omitted) . Therefore, if the partial action of a positive 
map on a state of a composite system spoils the positivity 
of the state, then the state must be entangled. Obviously, 
the Peres-Horodecki PPT criterion can be formulated as 
(f24|) . with A = T being the transposition operation. On 
the other hand, we note that realignment is not a positive 
map, and the related criterion involves the evaluation of 
the trace norm of the realigned state, which is in general 
not even Hermitian. 

One direction of the separability criterion based on 
positive maps can be applied in the space of matrices 
of moments to conclude that the starting state is entan- 
gled. Indeed, the reasoning at the base of the partial map 
criterion does not require any normalization and regards 
only the property of positivity. More explicitly: 

Criterion 6 Let A be a linear map preserving positivity 
of (infinite) matrices, and let M(p) be a separable matrix 
of moments, i.e., M{p) = ^ n p n M n (p A ) ® M n (p B ) with 
p n > 0. Then the (infinite) matrix resulting from the par- 
tial action of A, i.e., (id ® A) [M(p)] = J2 n Pn M n(p A ) ® 
A[M n {p B )], is also positive. 

Therefore, if we are given a matrix of moments M(p) 
for two modes and a positive map A and we find that 
(id® A) [M (p)] is not positive, then we conclude that the 
matrix of moments as well as the starting state are not 
separable. 

If there were a mapping between positive linear maps 
on states and positive linear maps on the corresponding 
matrices of moments, we could perhaps derive a general 
theorem of the Shchukin-Vogel type. Unfortunately such 
a connection, if existing at all, does not seem to be imme- 
diate. Transposition appears in this sense to be very spe- 
cial, since transposition of states translates simply into 
transposition of matrices of moments. Here, we will limit 
ourselves to the application of partial maps to truncated 
matrices of moments, so that we have the following: 

Criterion 7 //, for some f, there is a positive linear map 
A such that (id ® A)[M.-(p)] is not positive, then p is 
entangled. 

This Criterion is a direct consequence of the observa- 
tion at the basis of Proposition 1 and Criterion 6. Essen- 
tially, if one constructs a (sub)matrix of moments that 



preserves the separable structure of a state, and finds that 
the matrix of moments is entangled (using any arbitrary 
criterion, in this case linear maps), then one knows that 
the state was entangled. We remark that we are only able 
to establish a sufficient condition for entanglement (al- 
ternatively, a necessary condition for separability), con- 
trary to the analogous theorem for density matrices by 
Horodecki et al. [5J, which says that there always is a 
map able to detect the entanglement. 

We remark that in the case of transposition, which is 
defined for any dimension, the application of the map 
to the matrix of moments is equivalent to considering 
the matrix of moments of the partially transposed state. 
Therefore it is possible to directly focus on submatrices 
of the form Mf(p r ). On the other hand, in general, we 
may consider maps whose action is defined on finite di- 
mensions: consequently, we have to first take (properly 
chosen) submatrices Mj(p), and only then act partially 
on them to obtain M'~ = (id ® A)[Mj{p)}. This does not 
exclude that, after the action of the map, we may con- 
sider the positivity of an even smaller submatrix (M'~) T 
of the partially-transformed submatrix of moments. 



For example, one can apply non-decomposable 44 1 
maps to try to detect the entanglement of PPT entan- 
gled states. Classes of such maps were constructed for 
arbitrary finite dimension N > 3, e.g., by Kossakowski 
[29| , Ha l30ll. and recently by Yu and Liu [3l[ , Breuer (32| 
and Hall (331 ]. 

We are not able to provide examples of PPT bound 
entangled states, the entanglement of which is detected 
by applying positive maps on submatrices of moments, 
but the existence of such examples is not excluded. 

Furthermore, we stress that it may happen that a de- 
tection method based on an indecomposable map is able 
to detect more efficiently the entanglement of an NPT 
state than PT itself, e.g. it may be sufficient to consider 
smaller submatrices of moments. In any case, through 
the application of various indecomposable maps one can 
easily generate criteria for separability that are possibly 
independent from those obtained from PT. Indeed, as an 
important application of the proposed method we stress 
that it enables a simple derivation of interesting insepa- 
rability inequalities, e.g., 



2({N a N b ) + (NiN b )) < \(N a b) - (a?b)\ 



(25) 



which corresponds to the condition on the determinant 
of (f3"o]l obtained in the next subsection. 



A. Exemplary applications of positive maps 

The proposed method can be summarized as follows: 
First truncate the matrix of moments, i.e., M — > Mp 
then apply a positive map, i.e., — > M~, and check the 
positivity of the partially-transformed submatrix of mo- 
ments M'j. In turn, this amounts to considering positiv- 
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ity of submatrices (M'z) T , or, by virtue of Sylvester's cri- 
terion, to checking positivity of determinants det(M~) r . 
Thus, one can say that submatrices of partially trans- 
formed submatrices are considered. 

Here, we give a few examples of application of our in- 
separability criteria based on some specific classes of pos- 
itive maps applied to matrices of moments. 



the singlet state is given by i[3, —1, 1; — 1, 2, 1; 1, 1, 1] hav- 
ing negative determinant (equal to -1/4), which reveals 
the entanglement of the state. Analogously, the entan- 
glement of the partially entangled state \ip) = -i|(|00) + 

1 01) + 1 10)) can also be detected by (|2"9"j) . which is now 
reduced to (MUp)) T = §[4, -1,-1; -1,2,-1; -1,-1,1] 

with negative determinant (equal to -1/27). 



Kossakowski and Choi maps 



2. Breuer map 



The Kossakowski class of positive maps transforms ma- 
trices A = \Ag] iy x m in C N onto matrices in the same 
space as follows [29| 



A K [A] 



1 

TV" 



-TrA 



1 



N -1 



g-(Rx + nyTrA), (26) 



where '•' stands for the scalar product, k = -J (N — 1) /N , 
x = (xi)i, Xi = Tr {Agi}, and g = (g^i satisfying^ = g*, 
Tl {m9j} = S ij, Tr{.9i} = for i,j 



1. 



,N 2 



1. In 



our applications, we assume y = 0, R to be rotations 
R(9) e SO(N 2 - 1), and m to be generators of SU(iV). 
Note that the Ha maps [30| do not belong to (f2"6")l . In 
a special case for A = [A i A^ (3j the Kossakowski map is 
reduced to the Choi map [34j |. 

AchoifA] = -A + di&gQaAu + (3A 22 + ~/A 33 , 
7 An + aA 22 + f3A 33 , 
!3A ll + 1 A 22 +aA 33 ]) 1 (27) 

which is positive if and only if a > 1, a + /3 + 7 > 3 and 
l<a<2^>/?7>(2 — a) 2 , while decomposable if and 
only if a > 1 and 1 < a < 3 => f3j > (3 - a) 2 /4. We 
denote the resulting (unnormalized) matrix of moments 
shortly as 



M' f {p) = (id® A C hoi)[M f (p)] 



(28) 



It is worth noting that some bound entangled states can 
be detected [22j by applying to p the Stormer map [Hj], 
which is a special case of the Choi map for a = 2, j3 = 
0, 7 = 1 and of fl2B} for 9 = tt/3 and N = 3. 

Example. As an exemplary application of a positive 
map, let us apply the Stormer map to 9 x 9 matrix of 
moments M'-(p) for / = (1, a, a) ® (1,6,6). Note that 
the chosen map is indecomposable. For simplicity, we 
analyze only the submatrix (M~(p)) r for r = (2, 3, 7): 



Mn + M 22 -M 23 -M 27 
-M 32 M 22 + M 33 ~M 37 
-M 72 -M 73 M 77 + M 99 

l + (N a ) -(N a ) -(at 6) 

-(JV a ) 2(A a ) -(at 6) 

-(atft)' -(at6)* (NaN^ + im 



,(29) 



where Ma = (fjfj) are elements of the original (not- 



transformed) matrix of moments, Mz. Matrix 



for 



Our inseparability criterion for matrices of moments 
can also be based on the Breuer positive map defined in 
a space of even dimension d > 4 as follows [32[ : 



A B roucr[A] = lTlA-A-§[A], 



(30) 



where fl[A] = UA T U^ can be interpreted as a time re- 
versal transformation and is given by a skew-symmetric 
unitary matrix U . The latter can be constructed explic- 
itly as U = RDR T in terms of [33j]: 



D 



d/2-l 

e^ fc (|2fc)(2fc + l| - |2fc + l)(2fc|). (31) 



fc=0 



for any angles <pk and arbitrary orthogonal matrix R. 
Although antisymmetric unitary matrices exist only in 
even-dimensional spaces, the Breuer map can be gener- 
alized for arbitrary dimensions (see, e.g., [33|]). Thus, it is 
tempting to propose an analogous criterion by applying 
the Breuer map to a matrix of moments: 



M'f(p) = (id ® A BreU ex) [M f ~(p)} 



(32) 



and checking positivity of the transformed matrix M'~(p). 
It is worth noting that the Breuer map is a special case 
of the Yu-Liu positive map [31], thus even more pow- 
erful and computationally simple inseparability criteria 
for density matrices 3l|, 3^, H3] can also be applied for 
matrices of moments. 

Example 1. To reveal entanglement of the singlet state, 
let us first analyze a matrix Mi(p) of moments generated 

by some 16-element /. Antisymmetric unitary matrix 
U can, for example, be constructed as the anti-diagonal 
matrix 



U = 



o o o r 

o 10 

o-ioo 

-10 



(33) 



Then, by applying the corresponding Breuer map, one 
can easily get, from (|32p . the transformed 16 x 16 
matrix M'~(p) for arbitrary state p. This matrix re- 
veals, for example, entanglement of the singlet state 
for various choices of /, e.g.: /W = (l,a, N a ,a 2 ) <x> 
(l,6,A fc ,6 2 ), /CO = (l,a,A a ,l) <g> (l,6,A b ,l), or even 

;< 3 ) = (i,o, i,i)® (i,6, 1,1). 
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Note that and do not provide more informa- 
tion than (1, a, N a ) (gi (1, b, N b ) and (1, a) <g) (1, b), respec- 
tively. The matrices of moments corresponding to the 
former sets of operators contain redundant copies of the 
moments related to the latter sets, i.e., a repetition of 
an operator amounts to have a matrix of moments with 
repeated columns and rows. We considered such redun- 
dant sets of operators because Breuer criterion requires 
one of the subsystems to be at least 4-dimensional, but 
at the same time we wanted to emphasize that is possi- 
ble to detect (by means of Breuer's map) entanglement 
with fewer and fewer combinations of "independent" op- 
erators. We point out that provides for sure more 
information in general than f^ 2 \ and in turn the latter 
more than /( 3 ). 

The entanglement detection can be much simplified by 
analyzing the submatrix of M'~(fl) corresponding, e.g., to 

r=(2,5): 



Example 2. To reveal the entanglement of the Bell 
4,(100) + |11)), one can apply / = /« or 



state 



V2 



/( 2 ) and the Breuer map to be the same as in the for- 
mer example. Here, one can choose submatrix (M~(p)) r 
corresponding to r = (1, 6, 9), which reads as: 



M 2 ,2 - 


-M 3i3 


-Afi, e 


- M 3>8 


M2,10 - 


f M 3 ,ii 


-M e ,i 


-M 8 . 3 






-M 6)9 


- M 8 ,n 


M 10l2 - 


- M llj3 


-Af 9 , 6 


- M lli8 


Afio.io - 


f Mn,!! 



(38) 



For the analyzed Bell state, (J3SJ) yields det(Mj (1) (p)) r = 

det(M" j(2) (p)) r = —1/4 clearly demonstrating the en- 
tanglement. 

Thus, it is seen how new inseparability inequalities, 
corresponding to det(Mj(p)) r < 0, can be obtained by 
application of positive maps to matrices of moments. 



Mn 
-M 



Ma a -M 2S - M 



25 



M* 7 M 66 + M 



177 



VI. DETECTION OF BOUND 
(34) ENTANGLEMENT OF FINITE-DIMENSIONAL 
STATES THROUGH ANALYSIS OF MOMENTS 



where, as usual, My = (fjfj) are elements of the original 
matrix Ms-(p). For / = /W, matrix (f3~4"| reduces to 

(AfJ? (1) (p))r 

"(35) 

For the example of the singlet state, one gets 
(M| (1) (p)) r = [1, 1/2; 1/2,0], for which the determinant 
is —1/4. One can get even simpler criterion from (f3~4"|) by 
choosing / = f^: 



l + (a t2 a 2 ) -(a^b)-(a^ab) 
-(atft)* - (at 3 a&)* ((1 + N a )N a N b ) 



(Mj (2) (p)) r = 



2 (JV a 6) - (at 6) 

(JV a &t) - (abt) (iV a iV h ) + (7V a 2 AV 



.(36) 



Explicitly, for the singlet state, we have det(M~ (2) (p)) r = 
det[2, 1/2; 1/2,0] = -1/4. By contrast to /W and p\ 
matrix (|34|) for / = / < - 3 - ) is positive. Nevertheless entan- 
glement can be revealed by choosing a larger submatrix 
of M~ (3) (p) corresponding to r = (2, 5, 7, 8), which results 
in 



(M] (3) ( P ))r 







25* 2 J/1. 

y+ 2(N b ) y_ 
x* + y*_ z . 



(37) 



where jc± = ±(6) — (a^b), y± = ±(aN b ) — (Nb), and 
2 = ((JV a + l)Nb). For the singlet state, one again gets 
det(M| (3) (p)) r = -l/4. 

It is not surprising that one has to change submatrix 
(i.e. (|3"T|) instead of l[3"4"])). because for /( 3 ) less entries of 
the matrix Mf(p) contain independent information (ac- 
tually, only a 4x 4 matrix (corresponding to (1, a)®(l, b)) 
out of the larger 16 x 16 matrix (all the other entries are 
just repetitions)). 



The original SV criterion is based on partial transpo- 
sition, thus it cannot reveal PPT bound entanglement. 
On the other hand, it is known that the standard re- 
alignment criterion applied directly to the density matrix 
can detect entanglement of some bound entangled states 
[H HI El US VM- A question arises: Can PPT bound 
entanglement be detected by our realignment-based gen- 
eralized criterion? We have tested numerically some 
bound entangled states of dimensions 3x3 [IB, 36] , 2x4 
[H] , d x d H3, [H as well as infinite [13, GJ] , but we have 
not been able to detect entanglement by our generalized 
criterion. 

All numerical simulations suggest that the norms of 
reordered Mi satisfy the inequality v T ~ > uf or, equiva- 

lently, ||(M / -) r || > \\(Mj) R \\. If this observation is true in 
general, then the described realignment-based criterion is 
useless in detecting PPT bound entanglement. Neverthe- 
less, bound entanglement can be detected via moments 
with the help of the formula (see, e.g., [39j]): 



1 



{mi\p\m 2 ) = 



(-iy 



v / TOi!m 2 ! ^ J- 



E 



((a^) m2+j a mi+j ) 



which enables calculation of a given density matrix from 
moments of creation and annihilation operators. It is 
worth noting two properties: (i) The above sum is finite 
for finite-dimensional states (ii) Eq. (|39|) is not conver- 
gent for some states of the radiation field including ther- 
mal field with mean photon number > 1. The formula 
readily generalizes for two-mode fields as 



(mi,ni|p|m 2 ,n 2 ) 



OO 



((at) ma+J *a mi+J '(&t) n2+fc b ni+ *) 

(40) 
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Let us analyze a special case of (1401) for two qubits. 
Single-qubit annihilation operator is simply the Pauli op- 
erator given by a = o~ = [0, 1;0, 0], which implies that 
there are only four nonzero terms in sum (|40[) . We can ex- 
plicitly write two-qubit density in terms of the moments 
as follows: 



(NaN b ) 
(Nab), 

(aN b ), 
(ab), 



(N a N b ), 
(aN b ), 



(N a N b ) : 

(N a b), 



(at&t) 
(a^N b ) 
(iV a &t) 
{N a N b } 



, (41) 



coefficients c({ni,rrii}) that do not vanish. This corre- 
sponds to testing positivity of principal submatrices. 

For example, we show that Q implies the three-mode 
Hillery-Zubairy criterion [ll[ originally derived from the 
Cauchy-Schwarz inequality. By choosing / = (1, abc) 
(we use the notation introduced in Section |TTT| , one gets 
M f (p T ) = [1, (rfbc); (o&t c t), (N a N b N c )] , where N c = Sc 
and, analogously, N a and N b are the number operators. 
Imposing negativity of the determinant, one derives 



(N a N b N c ) < \(tfbc)\ 



(43) 



where N a = 1 - N a and N b = 1 - Nb. Matrix (J4TJ) 
can be partially transposed and realigned. All principal 
minors of p T are positive if and only if p is separable. 
The above simple example for 2x2 system was given to 
show the method only. To detect bound entanglement, 
one has to analyze at least 2 x 4 or 3 x 3 systems. For 
brevity, we will not present explicitly density matrices in 
terms of moments for these systems. Nevertheless, they 
can easily be constructed using ([40]) and then realigned, 
according to (fTUL to detect entanglement of some bound 
entangled states [2(J Hi], [22| ■ Finally, let us remark that 
there are drawbacks of the method: (i) it works if we 
know the dimension d < oo of a given state, (ii) Usually, 
it is simpler to directly reconstruct density matrix rather 
than to reconstruct it via moments. 



VII. A SIMPLE CONSTRUCTION OF 
MULTIMODE ENTANGLEMENT CRITERIA 

The two-mode SV criterion can readily be applied in 
the analysis of bipartite-entanglement of m-modes. For 
this purpose, one can define an m-mode normally-ordered 
operator 



f=f({ai}) = 



OO 



{ ni }=0 { mi }=0 



c({rii,mi 



})]>") 



f TTli 



(42) 

where for brevity we denote {n^} = {rii, n^, . . . , n m }, and 
similarly other expressions in curly brackets. As in the 
proof of Lemma[TJ we have that an operator X is positive 
semidefinite if and only if Tr {X ft f} > for every / as 
in (14"2"1) . To analyze how mode dj is entangled to all the 
other modes, it is enough to identify, in the reasoning 
followed in the previous sections, system A with the mode 
j and system B with all the other modes. Therefore we 
take a — aj, while normally-ordered powers tf ll b l2 are 
substituted by normally-ordered powers 

t(fci)i (fci) 2 J(k ] -i)i(k ] - 1 ) 2 



a, a 



tfe+i)i Akj+ih t(fc m )i (fc m ) a 
J + l j+1 ' ' ' m u m 



As in the two-mode setting, we may (and we will) analyze 
positivity of an operator X with respect to a restricted 
class of operators /, more specifically with only some 



which is the desired Hillery-Zubairy criterion [11] , i.e., 
a sufficient condition for the state to be entangled. By 
restricting the above case to two modes (say c = 1), one 
can choose / = (l,ab), which leads the Hillery-Zubairy 
two- mode entanglement condition , given by (|22p , as 
already shown in Q. By choosing a different function 
/, one can obtain new Hillery-Zubairy- type three-mode 
criteria. For example, let us choose / = (a, be) then 
M f (p r ) = [(N a ),(abc);(abc)*,(N b N c )}, which results in 
a sufficient condition for the three- mode entanglement: 



(N a )(N b N c ) < \ (abc)\' 



(44) 



In a special case, (j44|) is reduced to another two-mode en- 
tanglement condition of Hillery and Zubairy: (N a )(N b ) < 

S))| 2 , derived from the Cauchy-Schwarz inequality in 



VIII. CONCLUSIONS 

We have studied inseparability criteria for bipartite 
quantum states, which are given in terms of the matrices 
of observable moments of creation and annihilation op- 
erators, therefore generalizing the analysis by Shchukin 
and Vogel. Indeed, we have suggested (also by means 
of examples) that all the techniques originally developed 
to detect "directly" — that is, by considering the physical 
density matrix — the entanglement of states, can be de- 
ployed at the level of the matrices of moments. In doing 
this there are advantages — e.g., by considering an appro- 
priate submatrix of the matrix of moments one can apply 
techniques developed for finite dimensional system to de- 
tect the entanglement of infinite-dimensional systems — 
and disadvantages — e.g., while the separable structure of 
an entangled state is inherited by all properly constructed 
matrices of moments, it is not completely clear how the 
entanglement of the starting physical state gets encoded 
in the matrix of moments, and in some cases it may be 
difficult to choose the correct technique to detect it. 

In particular, we have proposed a new criterion based 
on realignment of elements of the moment matrices of 
special symmetry (i.e., corresponding to tensor product 
/s), as a generalization of the Rudolph-Chen- Wu realign- 
ment criterion applied for density matrices. Another re- 
ordering of elements of the moment matrices corresponds 
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to the partial transposition as in the original SV cri- 
terion. We have proposed another criterion based on 
positive maps applied to appropriate submatrices of mo- 
ments. We further observe that the formalism of matrices 
of moments can be certainly combined with the power- 
ful criterion invented in the finite-dimensional setting by 
Doherty et al. [40], in the attempt to detect, e.g., the 
entanglement of continuous-variable systems. How pow- 
erful this combination can be is nonetheless not evident 
or easily predictable, and we leave it as an interesting 
open problem. 

We have also discussed applications of the SV cri- 
teria to describe bipartite-entanglement of more than 
two modes. In particular, we have obtained the three- 
mode Hillery-Zubairy criteria originally derived from the 
Cauchy-Schwarz inequality, and derived new ones of the 
same type. 

As regards the confidence in the certification of entan- 
glement, if entanglement is verified within error bars for 
the matrix of moments (e.g., by considering the determi- 
nants of submatrices of the partially-transposed matrix 
of moments as in the original SV criterion), then entan- 
glement is certified for the physical state. This is true 
both in the case where error bars come from uncertain- 
ties in an experiment — from which the entries of the ma- 



trix of moments are obtained — or from numerical tools. 
We remark that here we are just considering certification 
of entanglement: in this paper we have not explored the 
relation between the degree of entanglement — as quanti- 
fied by some entanglement measure — of the physical state 
and the degree of entanglement of the matrix of moments. 

In conclusion, although it is an open question whether 
our criteria generalizing the Shchukin-Vogel idea are sen- 
sitive enough to detect bound entanglement, they enable 
to derive new classes of classical inequalities, which can 
be used for practical detection of quantum entanglement. 

Note added. After completion of the first version of 
our paper, the SV criterion was thoroughly applied to 
the multipartite CV case in [4l| . 
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